Phase stabilization of a frequency comb using multi-pulse quantum interferometry 
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From the interaction between a frequency comb and an atomic qubit, we derive quantum protocols 
for the determination of the carrier-envelope offset phase, using the qubit coherence as a reference, 
and without the need of frequency doubling or an octave spanning comb. Compared with a trivial 
interference protocol, the multi-pulse protocol results in a polynomial enhancement of the sensitivity 
0{N~^) with the number of laser pulses involved. We present specializations of the protocols 
using optical or hyperfine qubits, A-schemes and Raman transitions, and introduce methods where 
the reference is another phase-stable cw-laser or frequency comb. 



Quantum Physics has experienced a dramatic revolu- 
tion in the last decades, with unprecedented and uni- 
versally recognized [1] progress in the control and obser- 
vation of individual quantum systems. In this respect, 
trapped ions [2, 3] is one of the most mature setups, with 
unbeaten precision in the realization of single- [4] and 
two-qubit [5] unitaries and measurements [6, 7], closely 
followed by neutral atoms [8]. This spectacular progress 
underlies a number of "spin-off" , such as the characteri- 
zation of atomic properties using entanglement [9] or the 
development of quantum algorithms and protocols [10- 
12] for studying molecular ions. The synergy is even more 
advanced in the field of metrology, with accurate atomic 
clocks assisted by quantum gates [13, 14] or the use of 
atomic squeezing for enhanced magnetometry [15, 16]. 

Despite the exquisite precision of AMO systems, the 
times scales involved in control and detection (~ 10 /is to 
1 ms) have so far prevented exporting these techniques 
to study ultrafast processes. In this work, we show that 
the speed of AMO setups is not a limiting factor to accu- 
rately stabilize the carrier-envelope offset (CEO) phase 
of a frequency comb. CEO phase effects are relevant for 
few-cycle pulses, though effects in multicycle pulses have 
also been reported [17]. The first observation of such an 
effect was reported in the spatial asymmetry of above- 
threshold ionization (ATI) from Kr gas [18] and in the 
x-ray emission from Ne [19]. The direction of photocur- 
rents injected in semiconductors is also controlled by the 
CEO phase [20, 21] and more recently, the absolute CEO 
phase of single pulses has been measured in ATI experi- 
ments [22]. In general, the study of the CEO phase has 
been focused on its spectral components as these deter- 
mine the temporal evolution of the electric field within a 
pulse [23] , while only a few reports have addressed time- 
domain measurements of the relative phase of successive 
pulses in a train [24, 25]. The methods presented below 
follows this less-beaten path. 

Let us introduce the notion of "multi-pulse quantum 
interferometry" (MPQI), where an atom acts as a nonlin- 
ear, fast-response detector that efficiently measures the 
differences between ultrashort laser pulses. Modelling 
the atom-pulse interaction as a sequence of unitaries. 



{Ui}f^n through a suitable reordering of the pulses, ad- 
ditional gates and measurements, we build protocols that 
accurately determine the differences among the pulses, or 
the properties of the individual pulses themselves. Com- 
pared with cw laser interferometry, this approach prof- 
its from a polynomial enhancement of the sensitivity be- 
cause a single atom accumulates many interferometric 
events. 

A direct application of MPQI is the characterisation 
and stabilisation of a frequency comb. The comb [26, 27] 
is a device that produces an uninterrupted train of laser 
pulses with a fixed duration, r, and a regular spacing, T 
[cf. Fig. la]. Stabilising a comb is ensuring that the offset 
frequency, vq, remains a constant and well-known value, 
and that the spectrum is a collection of regularly spaced 
teeth with frequencies /„ = n/T -f vq. This problem was 
solved by Haensch and Hall [28, 29] in frequency space, 
interferometrically comparing different teeth in the limit 
of many pulses. Note that this requires a comb whose 
spectrum spans at least an octave, as some Ti:Sa lasers, 
or broadening the light with a nonlinear fiber. This sta- 
bilisation enables direct frequency comb spectroscopy, ac- 
curately revealing the (unknown) atomic level structure 
of neutral atoms [30, 31] and ions [32, 33]. 

We will rather work on the time-domain image of the 
pulse train. Here, the effect of the offset frequency, z/q, 
is to change the CEO phase from pulse to pulses, as in 
(t)n+i — 4>n — = vqT [cf. Fig. la]. To address the 
problem of comb stabilisation we will use MPQI, design- 
ing protocols that detect the phase difference between 
pulses with the greatest accuracy possible. The result- 
ing method does not require an octave-spanning comb, 
broadening or frequency doubling. It is thus useful for 
a wider variety of lasers, demands less power to perform 
the stabilisation, and it may profit from the ever-growing 
precision in atomic interferometry. 

Single-pulse unitary.- We start by determining the 
unitaries associated to each laser pulse and how they 
depend on the carrier-envelope phase, 0„. The inter- 
action of multi-level atoms with a frequency comb has 
been studied previously [34]. For simplicity and in order 
to produce analytical estimates, we will model this inter- 
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FIG. 1. (a) Electric field amplitude (solid line) and envelope 
(dashed) of a pulsed laser with period T and pulse-to-pulse 
phase difference Ac/i. (b) Associated spectra: a broad peak 
for one pulse (dashed) and a modulated comb for a pulse 
train with repetition rate /rop (solid). Note that the frequency 
offset uo depends on the pulse-to-pulse phase difference A0. 



action in the Rotating Wave Approximation (RWA) [35] 
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Here, m is the pulse index, s{t) > is the pulse enve- 
lope, id = 2Tri? is the comb carrier frequency, ujat is the 
atomic transition frequency {h — 1 throughout), there is 
an unknown (constant) phase (f>m for each pulse [36], and 
'^x,y,z are the Pauli matrices. The RWA works [37] for 
pulses which contain more than 30 periods of the carrier 
frequency, r > 30/z/, and allows us to explicitely write 
the pulse unitaries 

U,n = cos(0„O + i sin(0,„)a^„ = g-»0™<..^pe#™<T. ^ (2) 

in terms of the total Rabi flip angle of a single pulse, 6 = 
Jli^p s{t)dt, with a^^ = cos{(j)rn)<^x+sin{(j)m)o'y In what 
follows, we assume that the comb is almost resonant, 
uj — ujat, and has uniform intensity, that is 6^ — 0. 
These assumptions imply that we only need to stabilize 
the pulse-to-pulsc phase difference A(/). 

Multi-pulse unitaries.- We want a protocol that ef- 
ficiently detects the difference between the sequence of 
unequal pulses Utot = Y\!i=i and the ideal case . 
Let us first assume an ideal qubit, seeking an ordering 
of pulses with which the fidelity \ti{U^^Utot)\ decreases 
most rapidly with N. The simplest protocol (lA) applies 
N consecutive pulses [cf. Fig. 2a] with low intensity, 
6 ^ 1, on the qubit, which adiabatically follows the phase 
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Note how the pulse-to-pulse phase difference A0 de- 
creases the amplitude of the Rabi oscillations and can be 
measured. However, as we will show later on, the func- 
tional dependence on At/), implies a low sensitivity on the 
phase in practical implementations of the protocol. 

We can do much better by changing the intensity 
regime to = 7r/2 (protocol IB), where each comb pulse 
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FIG. 2. Different setups to measure the comb phase. A 
trapped atom interacts with either (a) one train of pulses 
or (b) two trains with a delay Td. Before and after the pulse 
train the atom is subject to one or more quantum gates and 
its state is interrogated at the end. (c) Experimental steps 
to perform a generalized Ramsey interferometer in order to 
estimate both A0 and 6. 



is capable of completely flipping the state of the atom. 
Under these conditions, for an even set of pulses we get 
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which for constant A(j) implies U^^^ = exp {—iNA(j)a^). 
Now A0 appears in a favorable position for the interfer- 
ometric detection with an enhancement proportional to 
the number of pulses, N. 

From Eq. (4) it is obvious that the sensitivity can be 
increased by maximizing the phase difference between 
consecutive pulses, (/)2fe — 4>2k-i- To take advantage of 
this, we design a set of protocols that split the original 
pulse train, extracting a sequence of N/2 pulses that will 
be delayed a time Td ^ T. This sequence of pulses will 
then be intercalated with the original one, as shown in 
Fig. 2b. The phases then read (/>2/c = kA(j) + A0 Ta/T 
and 4>2k-i = kA(j). Introducing this in Eq. (4) we obtain 
the unitary corresponding to protocol 2B: 



7-r(2B) 
'-^tot 



exp(-iCT^A0 X NTd/T) 



(5) 



which contains an enhancement factor, = T^/T. This 
is optimal with respect to any rearrangement of the 
pulses, using each pulse only once. 

Inter ferometry and sensitivity.- We must now trans- 
fer the information of the acquired phase to a measur- 
able quantity: the populations of the atomic states. For 
this, wc will complete the previous unitaries with ad- 
ditional operations and measurements that allows us to 
estimate A0 and 6. Out of 2M atoms, M are subject to 
the following steps [cf. Fig. 2c]: (i) initialize the atoms 
to the ground state, |0), (ii) apply a Hadamard gate, 
Uh ~ '^('^'■^ + ""z): (iii) apply a reference phase ^ onto 
the level (iv) let the atom interact with the comb 
as described before, (v) undo the Hadamard gate and 
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TABLE I. Experimental sensitivity, of a set of 2M two- 
level or three-level atoms to each of the protocols described 
in the main text. N is the number of pulses in a sequence, 
which in the delayed case are combined with A'' pulses from a 
later time, Td = NdT. 



measure the state of the atom, s e {0, 1}. The measure- 
ment outcome is described by the probability distribu- 
tion, Pi{s\9, A(j)). For the remaining AI atoms we skip 
(ii), obtaining the distribution, P2{s\9 , A(j)) . It is impor- 
tant to remark that we need no phase coherence between 
the comb and the lasers that implement the Hadamard 
gates, and that the reference phase, ^, is computed a 
priori to maximize the sensitivity of Pi^2 to the phase. 

The functions Pi and P2 convey all the information 
that is accessible in the lab: from the measurements of 
s in Pi and P2 experiments, one should compute dif- 
ferent estimators and use them to infer the values of 9 
and A0, with uncertainties ag and cta^- Using error 
propagation and the Fisher information we obtain fun- 
damental lower bounds and practical estimates [37] of the 
sensitivity (cr^J, and a^^) of each protocol (lA, IB, 2B, 
plus another combination, labelled 2A). As summarized 
in Table I, it is possible to build estimators of minimal 
variance for 9 and A0, which saturate the fundamental 
lower bounds. Moreover, we observe that all protocols 
but lA improve over the standard statistical sensitivity, 
-\/ M, thanks to the large number of pulses or to the use of 
pulses from well-separated times. In practical implemen- 
tations both N and N^i span several orders of magnitude, 
providing a sensitivity comparable to the state of the art. 

Three-level schemes.- In real atoms, if the qubit 
states and 1 are dipole-coupled by a comb, sponta- 
neous emission may severely limit the total duration of 
the experiment. One solution is to use dipole-forbidden 
transitions restricted in practice to the 9 <^ I regime. 
A more attractive configuration is the A-scheme shown 
in Fig. 3, where two long-lived states of the atom, |0, 1), 
talk via an intermediate level, |e). Applying combs or 
other lasers with orthogonal polarizations on the legs of 
the A-scheme, we can create effective Rabi oscillations 
between the qubit states, |0) and while still keeping 
the population of |e) so small that spontaneous emission 
is negligible. 

A simple way to minimize the spontaneous emission 
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FIG. 3. (color online) The m-th comb pulse interacts in a 
Raman setup with either (a) a cw laser signal or (b) another 
comb pulse. Controlling the polarization of the light and using 
the selection rules in atomic transitions we can ensure that 
each pulse or laser activates only one leg of the A scheme. 

is to turn the A- into a Raman scheme, detuning the 
lasers that couple |0, 1) with |e). Such Raman processes 
mix well with our algorithms. To start, if we have al- 
ready stabilized the phase of a cw laser, we can combine 
it with the pulses from the comb [cf. Fig. 3a]. This 
process enables an accurate determination of the pulse 
phase with respect to the cw source. The result will be 
a sequence of effective unitaries with some average Rabi 
angle, 9', and a pulse phase (f)'^ — (pm — iprci, where (/)rei 
is the phase of the stabilized source. The identifications 
9—^9' and 0™ — )■ 0J„ directly translate all protocols 
above to this new setup. Likewise, one may combine the 
frequency comb with a stabilized one [cf. Fig. 3b] and use 
our protocols to reconcile both sources. 

A more interesting use of Raman transitions is to 
achieve self-referencing with hyperfine qubits. For this, 
we repeat the scheme from Fig. 3b, combining two pulses 
from the same comb, but with a relative delay, T^. This 
amounts to a self-referenced interferometric scheme based 
on time shifts, not requiring frequency shifting nor shear- 
ing [23] . The phases of both pulses effectively combine in 
a nontrivial way in the unitary associated to the Raman 
process, 0^ = 4)^ - <j)m-Na = NdAcI). We can apply a 
sequence of iV pulse pairs with total angle which should 
optimally lay around N9' ~ 7r/4 

fj{lA, Raman) „ ^~iNdA4>cr^ ^iNO' ^iN^Acpa^ ^g-j 

and use Ramsey interferometry to measure both 9' and 
A0. A generalization of protocols 2A and 2B is also 
possible using a linear optics circuit with two delay lines, 
so that each atom is hit by pairs of pulses with alternating 
phases {4>m, 4>m-Ndi) and {4>m, 4>m-Nj2)- This leads to 
the sensitivities shown in the second half of Table I. 

Note that the use of Raman schemes demands the 
setup to be interferometrically stable up to a fraction 
of a wavelength. When a single pulse interacts with 
a two-level atom it does not matter whether the delay 
is a multiple of the comb period, or fails by a small 
amount, ST ^ - N^T {\5T\ < T). This is so be- 
cause only the phase of the pulse relative to the envelope 
enters the unitary and this only contains information on 
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v^NdT. However, in Raman schemes, where two pulses 
overlap in time, their relative delay appears as a new pa- 
rameter that influences the operations, not only in the 
accumulated Rabi angle, 6*', but also in the phase. In 
particular, the measured phase difference is actually [37] 
A(^' = A(/) + aj(5r, with a contribution due to the interfer- 
ometric path cJT, which must be separately stabilized. 

To remove the need for interferometric stability we can 
use a different A-scheme in which the comb only inter- 
acts with one leg, performing 7r/2 rotations. The un- 
perturbed and delayed pulses arrive closely in pairs, but 
without temporal overlap, implementing the sequence 
|1) ^ _e'("^'"-'^™-"d) Due to the lack of overlap, the 
delay errors drop and the effective operation is a phase 
gate in the qubit space. Spontaneous emission simply 
lowers the visibility and it is small because the excited 
state will only be populated a brief time Tf. = 0{t). De- 
noting by 7 the spontaneous emission rate of the |e) state, 
we may afford N ~ ~ log{e)/jTe pulses before we lower 
the visibility by e. Using a typical value of I/7 — 8ns, 
and a safe value — lOOps, spontaneous emission would 
be just 10% for about 200 pulses, sufficient to implement 
the last protocol in Table I. 

Errors.- An actual experiment will have sources of 
error which influence the atom's phase: ac Stark and 
Zeeman shifts induced by scattered light and stray mag- 
netic fields, and also small detunings of the laser. In all 
cases the effect is similar: a random term e{t)az makes 
the atomic levels fluctuate on a time scale which is much 
longer than r. Mathematically, this induces a time- 
dependent uncertainty in Aip of order x (tm+i — tm), 
where is the standard deviation of e(t) from its (zero) 
average, and tm the arrival time of each pulse. This er- 
ror can be cancelled using the spin-echo technique [38] 
or, more directly, in protocols 2A and 2B one can cali- 
brate the delays so that the pulses in a Raman pair arrive 
closely spaced, say 10 ps apart [39]. In this way, a pes- 
simistic ac Stark shift ~ 100 Hz induces an error of 
only 10^^ radians in the determination of A0. 

Another source of error is temperature: when atoms 
move between pulses, they sample the spatial variations 
of phase and intensity in the laser. We can eliminate 
both types of error [37] by (i) working in a Raman con- 
flguration which does not transfer net momentum to the 
atom and (ii) ensuring that the lasers are not tightly fo- 
cused. Both techniques also suppress any motion-qubit 
entanglement and the associated decoherence. 

The protocols discussed so far admit many differ- 
ent practical implementations: from neutral atoms to 
trapped ions. For concreteness, we will discuss this last 
setup, which has witnessed an astonishing progress in the 
fleld of ultrafast laser-ion interaction [40, 41]. The long 
coherence times in these setups, icoh ^ Is [42], allow in 
principle to consider trains of up to iV = t^oh/T ^ 10® 
pulses for a typical comb with /i-ep ~ 100 MHz [43]. 
Within the Raman schemes, using a single ion and one 



delay line, this already allows us to detect a pulse-to- 
pulse phase difference around ~ 10~® rad and thus 
calibrate the comb offset below vq < A(j)/T ^ 1/tcoh ^ 
1 Hz — a remarkable precision for 1 s interrogation time! 
The numbers improve, A(j) ~ 10~^^ using the same pulses 
in a A scheme, thus reaching a limit in which the pre- 
vious errors become relevant. The numbers decrease 
marginally, A^ ~ 10^ — 10^°, using shorter duty cycles 
~ 1 ms. These are optimistic sensitivities that go beyond 
what many existing lasers can hold: in existing combs, 
the phase may fluctuate violently after a shorter number 
of pulses. Our scheme cannot cope with that: we may 
only increase the delays Td and improve the accuracy up 
to a time scale at which these jumps happen. However, 
in practice such slow fluctuations set a limit on how well 
those combs may be stabilized by any method. 

Summing up, in this work we have presented an excit- 
ing line of research: the development of quantum algo- 
rithms for the characterization and stabilization of pulsed 
light sources using trapped atoms. Our work relies on 
the following key ideas. First, a single atom may ac- 
cumulate the effect of multiple laser pulses, computing 
their differences through the appropriate choice in the 
order of pulses, the intermediate gates and the final mea- 
surements. Second, there is a polynomial enhancement 
in the sensitivity of pulsed interferometric protocols with 
respect to conventional interferometry, arising from the 
combination of multiple events in a single atom. Third, 
these ideas can be used to detect changes in the phase of 
the pulses emitted by a frequency comb. This is possible 
because the unitary implemented by a single comb pulse 
is sensitive to both the intensity and the CEO phase, and 
not to the time at which the pulse arrives. 

The combination of all those ideas results in a variety of 
quantum protocols for stabilizing frequency combs. The 
schemes are particularly suitable for non-octave spanning 
combs with a low intrinsic phase noise, such as amplified 
high-power Ti:Sapphire laser systems, where significant 
phase noise is introduced by the amplification stages. In 
a near future, we expect to develop further pulse combi- 
nation and shaping protocols, and obtain formal results 
about the optimality of the algorithms, comparing their 
performance with existing linear and nonlinear interfer- 
ometric methods. Our protocols can be generalized also 
to work beyond RWA and to include and characterize 
other properties of the comb, such as pulse-to-pulse in- 
tensity fluctuations. We anticipate MPQI will enable new 
progress in flelds as diverse as ultrafast science, frequency 
metrology and direct frequency-comb spectroscopy, or co- 
herent control of molecular processes. 
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SUPPLEMENTARY MATERIAL 

In the following pages we provide more details for the 
following concepts which are spelled in the body of the 
paper: 

§A Justification of the Rotating Wave approximation 
and its influence in the following calculations. 

§B Study of the Raman transitions and three level sys- 
tems. 

§C Composition of unitaries for the different metrology 
protocols. 

§D Analysis of the sensitivity of the protocols using 
Fisher information and error propagation. 

§E Modelization of experimental errors and their in- 
fluence in the quantum gates. 

§F Practical implementation with trapped ions and es- 
timation of achievable sensitivies. 

§G Ultimate limits of these ideas based on pulse shap- 
ing techniques. 

§A. RWA AND CONVENTIONS 

We model the atom interaction with a single laser pulse 
in the semiclassical limit 

H ='^(T^ + s{t)cOs{u}t + (t),n)'J^. (7) 

Here, s{t) > is the pulse envelope, and d) = 2TrD is the 
comb carrier frequency. The evolution under this Hamil- 
tonian is described by a unitary operator that satisfies 
the Schrodinger equation i^wlt) = H{t)W{t). We are 
going to split out an evolution with the spin operator 
az as W{t) = exp{—iu]azt/2)U{t). This operator evolves 
now according to 

ij^U = ^^^^^^a, + .s(i)(e'^'+^^" +c.c.)(e'^V+ +H.C.). 

* (8) 

Assuming that the RWA is valid, in this Hamiltonian 
we can neglect the counter rotating terms, such as 
gz2wt-i-i0„^+^ and keep only those that are slowly varying. 
The result is i?RWA in the body of the paper. 

But when is the RWA valid? We have performed nu- 
merical simulations of the evolution of the qubit under 
Eq. (7), varying the duration of the pulse or number of 
oscillations it contains, as well as the intensity and de- 
tuning. In Fig. 4a we show the fidelity, F, of a resonant 
pulse, with a pulse area 9 = 7r/4, and a variable pulse 
length, T. The validity of RWA is also challenged by the 
inaccuracy of the control parameters, and in particular 
the driving frequency: as Fig. 4b shows, the unitary is 




10 20 30 40 50 60 
T(ji)/2n 




1.0 
0.8 

hi 0.6 
(Pm 0.4 



0.0 0.2 0.4 0.6 0.8 1.0 

FIG. 4. (a) Fidelity, F, between the unitary performed by the 
full model (7) and its description through the RWA (Huwa), 
as a function of the pulse duration, r, for a resonant pulse ui = 
ijJat- (b) Qubit excitation probability for detuned pulses with 
T — 30 X 2tt/ijJ (black) and 10 x 2n/u] (blue), and comparable 
Rabi frequencies. We show the exact solutions (solid) and the 
RWA (dashed), (c) Effective phase, (jjcS, of the unitary Um 
implemented by the full model (solid) and the RWA (dashed), 
as a function of the pulse phase, (j>m- 

affected by the detuning, and the differences between the 
RWA and the full model increase as the pulse length de- 
creases. In practice this is not a problem, for we expect 
the detuning of the comb to be smaller than 1%. As a 
final check, we show in Fig. 4c that the phase of the uni- 
tary is indeed proportional to the phase of the pulse. The 
main message is that for pulses above 30 oscillations, we 
are safe using the RWA Hamiltonian. 

§B. RAMAN TRANSITIONS 

Our Raman protocols are developed assuming that we 
can use ultrashort pulses to implement Raman transi- 
tions between two states, a and b, mediated by a third 
one, c, which remains unpopulated at the end of the 
pulse. (Due to the very short duration of the interroga- 
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FIG. 5. (a) Intermediate state population, Pc, and population 
of the state b, Pt, after a pulse of duration r, for cj — O.ScJat 
and Rabi frequency Q — O.OSSujat- (b) Evolution of the state 
populations during a pulse with Q = 0.076iUat, t ~ ^Qn/uj 
and uj = O.SoJat- (c) Relative phase between states b and a 
vs. the phase difference of the Raman pulse = (j>2 — in 
Eq. (9) for Q = 0.016tJ, a; = 0.98uJat and r = 4007r/a;. (Inset) 
Derivative confirming the monotonic behavior of (f>s. 



tion sequence, it is not necessary that c be unpopulated 
at all times, as is usual in STIRAP processes.) Impor- 
tantly, we need that such operations implement the same 
quantum gates and carry the same phase information as 
the original designs. We are going to discuss both re- 
quirements and how they are achieved. 

Note that Raman transitions with very short pulses 
have been demonstrated experimentally by the group of 
C. Monroe et al in a series of works that implement quan- 
tum gates with trapped ions and pulsed lasers [44, 45]. 



In these references, an interpretation based on Raman 
transitions induced by all the comb teeth is provided, 
but here we will discuss a different one. 

For us the key aspect of a Raman transition is the fact 
that the intermediate state, c, is completely depopulated 
at the end of the process. In order for this to happen, we 
need that the energy of the final state is similar to the en- 
ergy of the original one. Intuitively, this implies that the 
inverse of the duration of the process has to be smaller 
than 6 = LUat — i^, the detuning of the laser from the 
atomic transitions {a, b} o c, but larger than the differ- 
ence \u)ac — Wfccl- As shown in Figs. 3a-b, this qualitative 
appreciation remains true even for rather extreme cases. 
In those exaggerated plots, we see that pulses with a de- 
tuning d ^ 0.2iOat work fine even when they only contain 
10 — 20 oscillations of the laser. In this regime the excited 
state c is significantly populated during the pulse, but it 
has a population of less than 10~^ at the end. 

The other aspect we demand from these pulses is the 
fact that they must carry information on the phase of 
the laser. To check this, we analyze the interaction be- 
tween the three-level atom and the light using a simple 
Hamiltonian, 



H = s{t) cos{ujt + (pi) |c) (a| + H.c. 
-I- s{t) cos{ujt + (f)2) \c) {h\ + H.c 

+ UJat \c) (c| , 



(9) 



which under the RWA becomes 

HRWA{4>i,h) = s{t)e^' \c) {a\ + H.c. (10) 
^s{t)e"^^ \c) (&| -l-H.c 

+ {^at - W) |c) (c| . 

Note how HRWA{<t'i,4'2) is related to if_RVKA(0, 0) 
through a unitary transformation exp(— i(/)icr^^) in the 
{a, 6} subspace, with the relative phase 4>i = (1)2 — 4>i- I^i 
other words, according to the RWA the phase of the laser 
is mapped onto the relative phase between the states. 
The question is whether this behavior also follows from 
the original Eq. (9) . We have performed numerical simu- 
lations of the three level system in Eq. (9) and the conclu- 
sions are: (i) There is always a small deviation between 
the real phase and the RWA approximation, (ii) This 
deviation decreases with decreasing detuning, as in the 
two-level system, (iii) The actual phase experienced by 
the atom is a monotonic function of the laser phase, that 
is 4>s{4>i) grows with (j)i. These properties are exempli- 
fied in Fig. 5c for a case with 2% detuning, where the 
deviations from the RWA are small, below 1%, but the 
nonlinear behavior is clear in the inset. 

It would seem that, since we are striving for large accu- 
racies in the stabilization protocol, errors of 1% would be 
enough to discard the protocols. However, we have to re- 
member that we are not actually measuring the absolute 
phase, but the phase difference between pulses. Hence, 
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stabilizing i/is, which is a smooth, monotonic function of 
the laser phase, is equivalent to (and as accurate as) sta- 
bilizing 



§C. ANALYTICAL PULSE ESTIMATES 

We summarize some of the arguments in the body of 
the article regarding the composition of pulses. 



Pulse composition 

Let us start with the case of a sequence of pulses with 
^ <C 1 and a uniform carrier-envelope frequency mis- 
match, (j)^ = mAcj). The combination of pulses reads 

JV 

U^o^^ = n e'^''''e+''l''"''' (11) 

m=l 

JV 

= cos(6l) + isin(6») J2 e-^'t""''a^e^^'t"T'' 

m—1 

^ 1 + jge-'(^+^)^^ ''"^^,^,f a+ + H.c. + 0{6^) 
sm(A0) 

Note that in this context IS.(f) has the effect of a detuning 

and suppresses for long trains any excitation probability 
induced by the pulses. When we work around 9 = it/2 
we obtain instead 



(12) 



m=l 



If we assume that the number of pulses is even, we can 
use the anticommutation a^cTz = —(^z(^x and 

g-i<t>m<^z ^^Q+i<t>m''z g-i<t>m-icr^ ^^^+i(pm- 11^ z 

^ ^-2i{<l,^-4>^_,)a' ^ (13) 

recovering the formula 



Uiot = exp 
from the paper. 



N/2 

-2i^(02fc - ?f'2fe-i)o-' 



(14) 



Some optimality considerations 

We now prove that the sequence for protocol 2B (2 
sequence of pulses split from the original train with a 
time delay) is optimal when our only resource is the comb 
laser. As seen before, if we work around = Tr/2 we 
obtain the analytical formula 



Trim _ 
'^tot ~ ^xp 



-2i 



N/2 

E( 

fc=i 



l'2k 



(15) 



and our protocol accumulates phase quite fast, about 
0{NNd) for N pulses, where Nd depends on the delay. 
It is possible to prove that for any rearrangement of the 
same set of pulses (that is, with the same phases and 
intensity as before) this is the largest accumulation that 
can be detected. 

If (7 is a permutation for a certain arrangement of initial 
pulses, we can use the analytical expression for the arbi- 
trary product of a train of pulses with different phases to 
compute product of the unitaries after the permutation 



M 



<7(i) 



„2iE"i(-l)'''*'<A<.(i)<^ 



M 



(16) 



It is possible to find all permutations a such that 
they maximize \ J2iLi{~^y^^^'Pcr{i)\- Suppose the origi- 
nal pulses are ordered in terms of their carrier-envelope 
phase (pn. then it is quite straighforward to see how to 
construct optimal rearrangements of these pulses. Con- 
sider a a permutation of the first half of pulses and 
/3 a permutation of the second half of pulses. Then, 
the optimal set of rearrangements will be those formed 
by pulses labelled according to their carrier-envelope 
phase as {(j)a(i),(l)/3(N/2+i)}ie{i,...,N/2} or of the form 

{<p0{N/2+i),4'a{i)}iG{l,....N/2}- 

In particular, our proposed protocol corresponds to a 
and (3 being the identity permutation. This protocol ac- 
cumulates the largest possible amount of phase after the 
action of the pulses onto the ion. 



The fastest phase- accumulation protocol: phase referencing 

If we allow for more gates, performing unitaries 
in between the pulses, we can measure not only the 
phase difference, but also the total sum of the carrier- 
envelope phases J2i4'i- In order to do so, the new set 
of gates and unitaries, considered in order, would be 
{ax,U{(l)i), . . . ,ax,U{(l)N)}, for which the overall prod- 
uct is 



M 



M 



i=l i=l 
M M 

^J^^i't>i'^z^2^i<Pi<Jz ^J^^i't'i'^z^i't'iTz 



M 



where we use both GtC 



°J"^' a^. and al = Id. 



Note however this protocol demands gates in be- 
tween the pulses. Since the phase of these gates is stable, 
we can thus view this extra protocol as the referencing 
of the comb to the device that implements the cfx gates, 
which can itself be a laser or a microwave beam, in the 
case of hyperfine qubits. 
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§ D. FISHER INFORMATION AND SENSITIVITY 

We are interested in estimating the sensitivity of the 
interferometric protocols that we have developed with re- 
spect to changes in the parameters they depend on. A 
measure of the information that one can extract about 
one or several parameters from a given probability dis- 
tribution is the so-called Fisher Information [46, 47]. 

In our protocols, we want to estimate the intensity of 
each of the pulses 9 (considered constant throughout the 
whole experiment) and the pulse-to-pulse phase differ- 
ence A^. They will be related to some physical observ- 
ables which measure the population of the excited state 
after applying certain protocols. The precision of the 
parameters 6 and A0 is determined by the fluctuations 
of these observables and their variance can be obtained 
using standard error propagation theory. The Fisher In- 
formation will yield a measure of the available precision 
in the estimation of the parameters. Also, the variance 
of the estimation of a given parameter will be limited by 
the Cramer-Rao bound [48] , which sets the ultimate limit 
for the precision that we can achieve. 

Let us see how to compute both the Fisher Infor- 
mation and the Cramer-Rao bound. Let X be a sam- 
ple of observations with joint probability distribution 
given by P(X|k) depending on a vector parameter k = 
{ki,k2, . . ■ iki) and /i(k), a real valued function of k. 
Then, under suitable regularity conditions (see [46, 47]), 
for any unbiased estimator /i(X) of /i(k) 

Var(/i) > 6'^[I{k)]-^S 
where S is the vector of derivatives of /i(k), i.e. 




and the matrix /(k) is the Fisher Information matrix 
with (i, j)th element 

/„(k) = E(Si(X)Sj(X)) (18) 

where E(-) denotes the expectation value and 

5.(X) = ^logP(X|k). (19) 

In particular, if fc(X) is an unbiased estimator of the 
rth parameter kr, then, under the same regularity con- 
ditions 

Var(fcr) > Jrr(k), 

where JrrO<) is the rth diagonal element of the inverse 
matrix [/(k)]~^. 



We have both used standard error propagation theory 
and computed the Fisher Information as explained before 
in order to get fundamental lower bounds and practical 
estimates of each of the proposed protocols, which we 
have shown in Table I of the main text. According to the 
theory [48] , this implies the possiblity of finding estima- 
tors of minimal variance for both Acf) and 6. 

§E. EXPERIMENTAL ERRORS 

Dephasing 

We want to clarify with greater detail our estimates of 
the errors induced by small detunings and energy shifts. 
Our starting point is the Rabi model in the RWA, but 
now with an extra field, 

= ^a, + s{t) (e-^^'"<7+ + H.c.) . (20) 

e{t) which may be random but always so smoothly vary- 
ing, -^£{t) <C , that it may be assumed constant over 
a single pulse. If we make an interaction picture over this 
field, the result is that the phases of the pulses change in 
time 

= sit) (e-*^'"(T+ + H.c.) , (21) 

where the new phase depends on the arrival times of the 
pulses, t„i, as 

rtm+T/2 

0m - ^m-l — (Pm- 0m-l + / e{t)dt. (22) 

On average, such a field contributes to the phase differ- 
ence an amount which is of the order fj^{tm—tra-i), where 
CTg is the standard deviation of these random errors (as- 
suming they average to zero). We can thus upper bound 
the error in our phase estimation simply by decreasing 
the interval between pairs of consecutive pulses. This 
fine tuning is quite well suited for protocols 2A an 2B 
in the manuscript, where we are free to control the time 
separation between one pulse, t2k, and the pulse coming 
from the delayed beam, t2k+i- 

Temperature 

Temperature can also induce dephasing: if the atom is 
not still enough and it has time to move between pulses, 
it will see different phases of the pulse which depend on 
the distance traveled as 27rAa;/A. There are various ways 
to address this issue. We can make a simple estimate for 
a free atom in space, assuming that it is Doppler cooled. 
The temperature of the atom is 

keT ~ hV, (23) 
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where F is the natural linewidth of the cooling transition. 
Let us pessimistically assume that all this energy goes to 
the kinetic part, ^mv"^, giving us an average velocity 



2kBT 



2hr 



From this we can estimate the phase errors as 



A 



m — 1 



(25) 



Also pessimistically we will take F ~ 200 MHz, and a 
light atom such as Be, obtaining v ^ 5 m/s, which for 
a pulse separation of 10 ps gives 10"'^ (or actually a bit 
larger if we consider the additional velocity due to the 
photon recoil). 

The situation is very much improved when we arrange 
the lasers to work in a co-propagating Raman configura- 
tion such that there is no net momentum transfer to the 
atom. This is indeed a solution to the previous problem 
for, on each pulse, the spatially dependent phase which 
is acquired from one laser, kx, is cancelled by that of 
the other laser, —kx, making the whole process indepen- 
dent on the position of the atom. In this favorable cir- 
cumstance, the only effect that temperature may have is 
when the intensity of the laser varies spatially. However, 
by using a laser beam which is not too tightly focused 
and confining the atoms to a small region, the effect of 
this inhomogeneity may be safely neglected. 



§F. EXPERIMENTAL SENSITIVITIES WITH 
TRAPPED IONS 



As mentioned earlier, trapped atomic ions [2] consti- 
tute one of the most mature systems in the implemen- 
tation of Quantum Information Processing and Com- 
munication (QIPC) protocols and technologies. Preci- 
sion records have been achieved in the realization of 
single-qubit [4] and two-qubit [5] unitaries and mea- 
surements [6, 49], even reaching the threshold for fault- 
tolerant quantum error correction protocols [50]. In the 
last couple of years, fantastic progress in the controlled 
interaction between trapped ions and ultrafast lasers has 
been achieved by the group of C. Monroe at U. Maryland, 
with the demonstration of two qubit entanglement in the 
weak- field (0 <C 1), many-pulses regime [44], and also in 
the strong- field {6 ^ n), few-pulses regime [41], where 
logic gates faster than the trap oscillation period become 
accessible [40, 51]. Because of this, we think that the 
technology required to implement the phase-stabilization 
protocol that we propose is already available. 

Typical parameters for the frequency comb are listed 
in Table II. Pulses with a duration < 1 ps are nowadays 
easily accessible. On the other hand, one has to keep in 
mind that a pulse duration t effectively limits the possi- 
ble pulse delay times to > t in order to avoid an over- 
lap of the electric fields corresponding to different pulses: 



/.ep(MHz) 


r(ns) 


T (ps) 


Td (ps) 


100 


10 


10 


10-100 



(24) TABLE II. Typical parameters of a frequency comb [40, 41] 



depending on their polarizations, this may lead to several 
unwanted effects, from excitation of motional sidebands 
to total cancellation of the Raman transition [40]; in ei- 
ther case, the action of a pulse pair on the qubit would 
still be described by a unitary transformation Ui, but not 
the ones we have written down earlier, so that our model 
would break down. Therefore, we will stick to a comb 
with pulses of 1-10 ps. 

In the following, we present details of our calculations 
to estimate the achievable sensitivity enhancements for 
the protocols introduced in the main text of the article. 
At the end, we present an estimation of the ultimate 
sensitivity limit that can be reached with pulse shaping 
techniques. 



Sensitivity enhancements with two-level protocols 

Protocols lA to 2B consider direct transitions induced 
only by the comb laser that we want to study. In prac- 
tice, there are two ways that this can be achieved: dipole 
and quadrupole transitions. Dipole transitions are, for 
instance, the ^5'i/2-^-Pi/2,3/2 lines in Yb+ [45]. These 
transitions have a typical linewidth of a few tens of MHz, 
which is comparable to /rop- This implies that the time 
between consecutive pulses could be shorter than the life- 
time of the excited state of the ion. A possible solution 
to this problem will be presented later on in §G. 

An alternative is to rely on quadrupole transitions, as 
provided by the Ca"*" ion using as qubit states the elec- 
tronic states IS'1/2) and I-D5/2) [3]. The excited level now 
has a radiative lifetime Tj-ad ^ 1 s which is favorable to 
implement our ideas. The downside of quadrupole tran- 
sitions is their lower coupling strength, which demands a 
more powerful laser to excite them. In practice, depend- 
ing on the laser, this might imply that we have to work 
in the limit <C 1 (protocols lA and 2A) but we will 
forget this in following discussion. 

Following Table II, let us consider a frequency comb 
composed of 10-ps pulses with /rep = 100 MHz. This 
pulse duration is much shorter than the trap oscillation 
period (1 fj,s for a typical Wtrap = 27r x 1 MHz rf Paul 
trap) and allows us to disregard the motional state of 
the ion in the trap [cf. Eq. (7)] as well as its micromo- 
tion, which may affect the performance of coherent pro- 
tocols at longer times [45]. Such a frequency comb, and a 
typical ion coherence time for the electronic qubit states 
of '*°Ca+, Tcoh ~ 10 ms [3], allow up to Tcoh/rep ~ 10^ 
pulses to go through the ion before decoherence becomes 
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relevant. Using protocol IB, this leads to an enhance- 
ment of the sensitivity by a factor xiB = N = 10^, 
which translates in a stabilization of the frequency off- 
set down to Svo — f-ccp/N ^ 1/Tcoh ~ IkHz. (In ad- 
dition, a probe time of 10 ms combines well with ion 
state-detection times ~ 10ms for an optimal duty cycle.) 

This result can be improved by applying the pro- 
tocols with two pulse sequences (2A,2B). To be spe- 
cific, we can pair iV = 5 x 10^ pulses with delays 
Nd = 5 X 10^ and reach, with protocol 2B, a resolution 
'5i^o = /rcp/(iV7Vd)=0.4mHz. 

Again, we remark that the numerical estimates in the 
main text of the paper and this supplementary mate- 
rial take into consideration only the coherence properties 
of trapped ions for the stabilization an "ideal frequency 
comb" , and technical issues inherent to currently avail- 
able combs are not included in the calculations. 

Sensitivity enhancement with self-referenced Raman schemes 

Use of a Raman scheme lifts the restrictions related 
to the excited-state lifetime of the qubit as spontaneous 
decay is of no concern. Such a scheme has been imple- 
mented with various systems, e.g., Yb^ with a qubit de- 
fined by the ^5*1/2 hyperfine states \F — 1, nip = 0) = |1) 
and \F — 0,mp = 0) = |0), which are split by ujat — 
12.642815 GHz [40, 41]. For these states, coherence times 
larger than 1 second have been measured [42]. 

Let us consider a pulse train of 1ms, which provides 
10^ pulses at f^^p = 100 MHz, and let us split this train 
on two lines. We seek to maximize the phase difference 
between them. To this end, we consider the available 
10^ pulses into sets of lO'' and keep the first set. Si = 
{1, . . . , lO''}, the set 52 = {lO"* + 1, . . . , 2 x 10^}, and the 
last set, ^3 = {9 X 10^ + 1, . . . , 10^}. The first set, Si, 
will be further split in two, so that half of the pulses are 
paired with those in S2 {Ndi = 10*) and the other half 
with S3 {Nd2 = 9 X lO''). Then, this optical setup yields 
a sensitivity enhancement of order x = N\Nd2 — Ndi\ = 
8 X 10® or Si^o ^ 0.1 Hz. If we allow ourselves a longer 
interrogation time of 1 s, the figures would improve down 
to an amazing precision of 10~^ Hz. 

We note that these high sensitivies are achievable al- 
most independently of the underlying physical system 
used for the qubit: taking into account the continuum 
spectrum of each pulse, the only requirement is the prox- 
imity of Wat and (D, a feature that can be engineered, and 
coherence times which are experimentally available. 

Recursive refinement for large pulse-to-pulse phase shifts 

In our studies we have found that the sensitivity of 
our metrology protocols can be written in the form 
^ x{N)\/M, where the enhancement factor x(^) 



arises from a clever accumulation of the phase. In prac- 
tice, for a non-stabilized frequency comb with a large 
and an excessive number of pulses, the total accumulated 
phase, x(iV)A(/), will wrap around the maximum measur- 
able value, TT, precluding a unique determination of Aip. 

The appropriate way to deal with this situation is to 
do an iterative refinement of the phase measurement. 
As an example, let us consider a fiber-based frequency 
comb: these devices have an intrinsic width of the off- 
set frequency of about 200 kHz. This means that, for 
/icp = 100 MHz, the phase ipm may wrap around tt in 
about N = 500 pulses. Hence, on the first iteration, it is 
meaningless to interrogate the laser for much longer than 
a few fis. This iteration allows us already to achieve a 
precision in A(f> of order •\/M/500 with protocol IB. This 
initial value can be used to fed back to the laser setup, to 
lower A(j) and, on the next iteration, use a larger number 
of pulses. 

Continuing with this example, a similar iterative re- 
finement using protocol 2B and a fixed interrogation time 
~ 1 /iS, would lead to an accuracy in the comb offset fre- 
quency of SiyQ = f rep/ (NNd) = /rep/2502 = 3Hz using 
only one ion. 



§G. ULTIMATE PRECISION LIMITS WITH 
ADVANCED PULSE-SHAPING TECHNIQUES 

The protocols discussed so far achieve a great efficiency 
thanks to the number of pulses in a given interrogation 
time and possible delays among them. Note however, 
that the comb is mostly "empty": between every two 
pulses of about 10 ps, there is a waiting time ^ 10 ns 
in which the ion is idle. It would seem that this empty 
time, combined with the coherence rates of the ions, set 
the ultimate limits for precision in our setup. However, 
if the laser has enough power, we can engineer a clever 
scheme to fill these empty gaps, increasing the effective 
repetition rate of the ion-laser interaction. 

The trick here will be to "compress" the pulses so that 
a minimal time elapses between the end of one pulse and 
the beginning of the following one, but without modifying 
the phase of any one pulse. Such "compression" could be 
realized with an optical setup as depicted in Fig. 6b. The 
key ingredient in this setup is an optical device which we 
call Beam Splitter and Delayer (BSD) that, given an in- 
tense ultrashort pulse, extracts a train of n replica pulses 
separated by a very short time At [cf. Fig. 6a]. (An al- 
ternative BSD optical setup producing 8 replicas of an 
initial pulse has been recently implemented in [41].) 

We will discuss these ideas in a particular application: 
doing metrology of the comb with a dipole transition. In 
this case the qubit of choice will not satisfy the condition 
Tcoh ^ T. Consider for example Ca''" ions using the 
dipole-coupled 1 6*1/2) and l-pL/2) states for which Trad ~ 
7 ns [52] . In this setup we can still reach high precisions 
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FIG. 6. (a) A beam splitter and delayer (BSD) can be implemented by shining the laser pulse train into a semitransparent 
mirror facing a perfect mirror (shaded area). The incident angle of the beam determines the number, n, and delay, At, between 
the outgoing replica pulses, which are focused by a lens (L) onto the ion. (b) Optical setup with a delay line (dashed box) and 
two BSDs with delay At = 4r. This transforms the incident pulse train into a new train with interleaved pulses with a fixed 
phase difference between them and a total duration 2(n — l)r kT. 



taking a relatively long pulse train of duration ^ Tcoh as 
long as we ensure that all the pulses pass through the ion 
within a short time ^ Tcoh- On the other hand, wc must 
still fulfill the requirement that different pulses do not 
overlap in time, that is At > 2r, which restricts us to use 
sets of up to n < min{Ti-ad, Tcoh}/(2T) ^ 7 ns/20 ps — 350 
replica pulses. To be concrete, let us use the setup in 
Fig. 6b to pick up two pulses with a relative delay of 
Td = NdT ~ lOfis — this corresponds to pulses 1 and 
k = Nd = 10^ in the previous figure. The pulses will go 
through the BSD and be recombined, alternating replica 
pulses from each line. For a conservative n — A (not 
to lose too much power in each replica), At = 4t, and 
Nd « 1000, we obtain a phase sensitivity enhancement 



by a factor nNd « 4 x 10^. 

The same ideas can be applied to the Raman scheme 
by ensuring that the replica pulses from both lines arrive 
simultaneously to the ion. The result is an enhancement 
of the sensitivity by an additional factor n on top of the 
formulae derived in the main text of the paper. 

We finally note that the very short probe times con- 
sidered here, allow for the recollection of a large set of 
statistical data in a very short time. Together with the 
large sensitivity enhancements calculated, the presented 
schemes appear as very competitive protocols to mea- 
sure and stabilize the carrier-envelope offset phase of fre- 
quency combs without the need for octave-spanning spec- 
tra. 



